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ABSTRACT 

Linear analysis of the formation of protostellar cores in planar magnetic interstellar clouds shows 
that molecular clouds exhibit a preferred length scale for collapse that depends on the mass-to-flux 
ratio and neutral-ion collision time within the cloud. We extend this linear analysis to the context 
of clustered star formation. By combining the results of the linear analysis with a realistic ionization 
profile for the cloud, we find that a molecular cloud may evolve through two fragmentation events in 
the evolution toward the formation of stars. Our model suggests that the initial fragmentation into 
clumps occurs for a transcritical cloud on parsec scales while the second fragmentation can occur for 
transcritical and supercritical cores on subparsec scales. Comparison of our results with several star 
forming regions (Perseus, Taurus, Pipe Nebula) shows support for a two-stage fragmentation model. 
Subject headings: diffusion — ISM: clouds — ISM: individual (Perseus, Pipe Nebula, Taurus) — ISM: 
magnetic fields — magnetohydrodynamics — stars: formation 



1. INTRODUCTION 

Molecular clouds are observed to have complex mor- 
phology existing on a wide range of sca les (e.g., 
iMcKee fc Ostrikenl2007tlGoldsmith et al.ll2008D . and ex- 
hibit numerous clumpy and filamentary structures. They 
are the birthplaces of stars, and the observed distribu- 
tion of young stellar objects (YSOs) also yields a wide 
range of object densities. The YSOs are measured to 
be anywhere from lo osely clustered to highly clustered 
(jMegeath et al.ll2009fl . 

Many molecular cloud maps do however reveal an in- 
teresting pattern of dark clouds within which parsec- 
scale clusters are being formed. The dark clouds are 
often themselves separated by several pc. Even the 
Taurus molecular cloud, often associated with so-called 
"isolate d" star formation, demonstrates several dark 
clouds ( Onishi et al. 1998) of pc-scale and masses ~ 
100 Mq that are separated by several pc, each con- 
taining embedded weak clusters of YSOs. It was sug- 
gested long ago that the dark clouds within a molec- 
ular cloud complex m ay be the result o f some early 
fragmentation process (iDutrev et aT]|199H iGaida et al.l 
11984 iSchneider fc Elmegreenl I1979D . and subsequently 
the suggestion was made that the fragmentation within 
the dar k clouds may be controlled by the local Jeans 
length (jHartmannl I2002T ) . Based on the observed mor- 
phology of clustering within parsec-scale clouds, there is 
indeed a case to be made for a two-stage fragmentation 
process. The fragmentation within the dark clouds also 
shows a possible scaling with Jeans length, e.g., the typi- 
cal separation between cores in the Taurus dark clouds is 
0.25 pc, consistent with the Jeans length corresponding 
to a column density associated with the visual extinction 
A v ~ 5. 

The question then arises as to what controls the frag- 
mentation process at the parsec scale, leading to the 
formation of individual dark clouds within a molecu- 

nwityk@uwo.ca (NDB); basu@ uwo.ca (SB)| 



lar cloud complex. Since these regions have nonthermal 
line widths, the mechanism of turbulent fragmentation, 
in whi ch turbulence plays a dominant role, is often in- 
voked (jPadoan et al.lll997t IGammie et all 120031; IKlessenl 
120011 ). However, turbulent fragmentation models make 
no unique prediction for any observed molecular cloud. 
The usual technique of running simulations with power 
spectra of the form P(k) oc k~ n , with n > 0, means that 
much of the power is contained in the largest scale of 
the simulation, leading to the largest structures forming 
on the chosen size scale of the simulation region. In this 
paper, we look for an explanation for large initial frag- 
mentation scales that does not depend on the scale of 
the modeled system. We employ a linear analysis that 
assumes a uniform background state, so it does suffer 
from the standard "Jeans swindle" that is well known in 
the nonmagnetic analysis. We apply the analysis results 
using thermal line widths, without using the ansatz of 
replacing the isothermal sound speed with a nonthermal 
line width. We are in effect assuming that the fragmenta- 
tion of a diffuse cloud is occurring in the dense midplane 
that has thermal line widths, while an extended diffuse 
envelope can have nonthermal line widths. Such a pic- 
ture is consistent with the modeling of th e propagation of 
turbulent m otions in a stratified cloud by Ku doh fc Basul 

<l200af2006l)- 

iCiolek fc Basul (|2006) carried out linear analysis of 
gravitational instability in a flattened layer that is par- 
tially ionized and threaded by a magnetic field that is 
initially perpend icular to the s heet. This is based on 
earlier work by iMortonl (jT991) . The results show a 
dramatic dependence of the preferred scale of fragmen- 
tation upon the ambient mass-to-flux ratio, if this ra- 
tio is close to the critical value. The time scale for 
growth of the instability also exhibits a dramatic change 
around the critical mass-to-flux ratio, with a step- like 
drop in value as the mass-to-flux ratio goes from slightly 
subcritical to slightly supercritical. The time scale of 
growth of gravitational instability in the subcritical case 
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is commonly referred to in the literature as the "am- 
bipolar diffusion time" . A linear analysis of an infi- 
nite uniform medium with partial ionization and allow- 
ing for ambip olar diffusion has also been performed by 



Langcr (1978) and recen t ly more comprehensively by 
Mouschovias et al.1 ()2011h . iKunz fe Mouschoviasl (2009) 



have also used this linear analysis to generate a core 
mass function to fit the one observed in Orion. In 
this paper, we inco r porat e the new theoretical results 
from lCiolek fcBas u (2006) regarding the fragmentation 
of partially ionized and magnetized gas, focusing on the 
fragmentation of a flattened layer, as might be an ex- 
pected early configuration if the magnetic field strength 
is significant, or if a cloud is formed by large scale flows. 
Furthermore, structured objects like sheets or filaments 
exhibit a preferred scale of fragmentation, unlike uni- 
form media, for which the largest scale always has the 
greatest growth rate. In a sheet-like situation, initial 
small-amplitude white noise perturbations result in the 
preferred scale (possessing the minimum growth time in a 
linear analysis) e ventually growing to dominate the non- 
linear evolution (|Basu "eTaLll2009H) . 

Another dramatic variation of the ambipolar diffu- 
sion time scale in a molecular cloud will occur be- 
cause of the strong dependence of ioniz ation fraction 
on column density of a molecular cloud (HVIcKed IT989t 
ICiolek fc Mouschoviasll99llRuffle et al.lll998h . The ion- 
ization fraction can drop in a step-like manner as the 
column density of a cloud increases beyond a threshold 
value. The long ambipolar diffusion time scale in the 
UV-ionized region has been invoked as the re ason for 
the r elatively low Galactic star formation rate (|McKed 
1989), and the transition to lower ionization due to cos- 
mic rays has been invoked as a reaso n for lower turbu- 
lence in high column de nsity regions (jRuffle et al.llT99i 
iMyers fc Lazariadfl998l ). The effect of this drop m ion- 
ization upon the fragmentation process itself has not pre- 
viously been investigated. 

Based on the above findings, one can envision two 
dramatic events during the assemblage of a molecular 
cloud. If we assume that a molecular cloud is assem- 
bled fi^m_Fri_claud_inaterial that is generally subcriti- 
cal (jHeiles fc Trolandll2005l) and ionized by UV starlight, 
then a cloud may resist fragmentation (due to the very 
long ambipolar diffusion time) until flows primarily along 
the magnetic field lines raise the mass-to-flux ratio to 
slightly above the critical value. At this point, the trans- 
critical (but still UV ionized) cloud may form large frag- 
ments as predicted by the linear theory. As these frag- 
ments develop, they will also become more supercritical, 
since their evolution is partially driven by neutral-ion 
drift. Once the column density in the fragments also 
crosses the column density threshold for the transition 
to cosmic ray dominated ionization, the fragmentation 
length and time scales will drop dramatically, and a sec- 
ond stage of fragmentation may be possible. This is the 
scenario that we investigate in this paper, by combining 
results of linear analysis and ionization balance. 

In Sect ion 2 we outline our model and revi sit the linear 
analysis (|Mortonlll99lt ICiolek fc Basull2006l ) with a spe- 
cific focus on our chosen model parameters. In Section 
3, we use the linear analysis results to discuss the feasi- 
bility of core formation by subclumps forming within a 
larger contracting diffuse clump. We do this by compar- 



ing numbers from the model with core and clump length 
scales from observed star-forming regions. A further dis- 
cussion is given in Section 4 and a summary is presented 
in Section 5. 

2. PHYSICAL MODEL 

We consider clump/core collapse within partially ion- 
ized, isothermal, magnetic interstellar molecular clouds. 
These clouds are modeled as a planar sheet with infi- 
nite extent in the x- and y-directions and a local verti- 
cal half-thickness Z. The nonaxisymetric equations and 
formulations of our m odel have been described in de- 
tail in several papers (|Ciolek fc Basul[2006l : iBasu et al.1 
l2009alf bT). The volume density (p n ,o) is calculated from 
the vertical pressure balance equation 



PnfiC f _ 



Pc 



(1) 



where P cx t is the external pressure on the sheet, a n fl 
is the initial uniform column density of the sheet and 
c s = {kgT /uin) 1 / 2 is the isothermal sound speed; ks 
is the Boltzmann constant, T is the temperature in 
Kelvins and m n is the mean mass of a neutral particle 
(m n = 2.33 amu). The evolution equations also include 
the effect of ambipolar diffusion, which is a measure of 
the coupling of the neutrals with the ions, and by exten- 
sion the magnetic field. This is quantified by the time 
scale for collisions between neutrals and ions: 
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Here, rrn is the ion mass, and (aw)iH 2 is the neutral-ion 
collision rate. The typical atomic and molecular species 
within a molecular cloud are singly ionized Na, Mg and 
HCO, which have a mass of 25 amu. Assuming collisions 
between H2 and HCO + , the neutral- ion collisio n rate is 
1.69 x 10~ 9 cm 3 s" 1 (|McDaniel fc Masonl [l97l . These 
collisions transfer information about the magnetic field 
to the neutral particles via the ions that are bound to 
the field lines. The threshold for whether a region of a 
molecular cloud is stable or unstable to collapse is given 
by the normalized mass-to-flux ratio of the background 
reference state, 



Mo 



(3) 



where (2-7TG 1 / 2 ) 1 is the critical mass-to-flux ra- 
tio for gravitational collapse in the adopted model 
(jNakano fc Nakamural H978t ICiolek fc BasuT 120061) and 
-B re f is the magnetic field strength of the background ref- 
erence state. Regions with fiQ < 1 arc defined as sub- 
critical, regions with fio > 1 are defined to be supercrit- 
ical, and regions with /i ~ 1 are transcritical. In the 
limit of flux freezing, the time scale for collisions (r„i) 
is zero, implying frequent collisions between the neutral 
particles and ions, which therefore couples the neutrals 
to the magnetic field. Under these conditions, subcriti- 
cal regions are supported by the magnetic field and only 
supercritical regions can collapse. For nonzero values 
of r n i, a finite time scale for collisions between neutral 
particles and ions exists which is inversely dependent on 
the ion density. The variation in ion density within a 
cloud is linked to the type of ionizing radiation available 
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at different depths with i n the cloud. This is shown in 
Figure 1 of iRuffle et al.l (|1998| ). At low extinction, i.e., 
on the edge of the cloud, the ionization fraction is large, 
and declines steeply beyond a threshold column density 
due to shielding of the inner regions from the UV ra- 
diation that ionizes the outer layers. Ionization within 
these inner regions is predominantly due to cosmic rays. 
Ambipolar diffusion, i.e., neutral- ion slip, can lead to 
redistribution of mass relative to magnetic flux within 
a molecular cloud and cause gravitationahy unstable re- 
gions to develop within subcritical regions. These regions 
then have the chance to collapse. 

The model we use is characterized by several di- 
mensionless free parameters. These include a dimen- 
sionless form of the initial neutral-ion collision time 
{Tnifl/to = 2irGa ni oT n i.o/c s ) and a dimensionless ex- 
ternal pressure (P cx t = 2P cxt /7rGof i ), in addition to 
the dimensionless mass-to-fiux ratio /io that was intro- 
duced earlier. We normalize column densities by <7 n ,0j 
length scales by Lq — c 2 /2nGa n! o and time scales by 
to = c s /2irGa n fi. Based on these parameters, typical 
values of the units used and other derived quantities are 
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where n n ^ is the initial neutral number density and Xifl 
the ionization fraction. Note that T n i t o is inversely pro- 
portional to Xi,o f° r a fixed cr nj o- For our analysis, we 
assume a dimensionless external pressure P ex t = 0.1 and 
a temperature T = 10 K. 

2.1. Linear Analysis 

Assuming the unperturbed zero-order or background 
state of the model is static and uniform, th e magneto- 
hydrod ynamic equations for a thin sheet (see iBasu et al.l 
2009a) can be linearized to first order for any physical 
quantity via 

f(x,y,t) = f + 6f 

= fo + 5f a e t{k * x+k v y ~" t \ (9) 

where /o is the unperturbed background state, Sf is the 
perturbation, Sf a is the amplitude of the perturbation, 



kx 5 ky ; cind k are the x-, y-, and z- wavenumbers respec- 
tively such that k 2 = k 2 + k 2 and lu is the complex an- 
gular frequency. The perturbations are necessarily small 
such that \5f a \ <C /o- With this assumed perturbation, 
d/dt — s> — iu), d/dx —> ik x , and d/dy —> ik y . Apply- 
ing the perturbation to the dimens ional equat i ons for a 
model cloud (Equations (2) - (4) of IBasu et all (|2009al) ). 
and retaining only the first order te rms, one obtains, in 
the manner of iCiolek fc Basul (|2006| ). 

+ k y C S SV nt y, (10) 

k 

+ y {2nGa nfi ^ 1 + kVl Q ^)8B z ^ (11) 
k 

uJc s Sv n . y = -^(C 2 S0 k - 2itG(j n fi)8o nfi 

+ ^(27rG<7„,o^o 1 + *V2 i0 /io)*B»,eq, (12) 
ujSB z eq = k x c s 8v n x H k v c s Sv n y 

-iT nifi (2irG<j nfi fj,Q 2 k + k 2 V^ )SB ZyCq , 

(13) 

where v UiX and v n . y are the velocities of the neutral par- 
ticles in the x- and y- directions respectively, P 2j0 q is 
the magnetic field strength at the midplane of the sheet, 
G^o is the local effective sound speed, where 



n 2 2 [3Pext + (7T/2)G< ] 2 



'n,0 



[Pext + (ir/2)Ga 2 



2 12 s' 



(14) 



Va,o is the Alfven speed which is related to the mass-to- 
flux ratio (/io) v i a 



vi, = 



B: 



= 27rG(T„ j oA < o 2 Zq, 



(15) 



and Zq is the initial half-thickness of the sheet. 

A mode is unstable if the imaginary part of the com- 
plex frequency wim > 0. The growth time of such an in- 
stability is T g = I/cjim- The dispersion relation is found 
to be 

(oj + i6){uj 2 - C 2 Sfi k 2 + 2irGa n , k) 

= uu{2nG<T nfl k^ 2 + k 2 Vl Q ), (16) 

where 

9 = T mfi {2^Ga nfl k^ 2 + k 2 Vi ). (17) 

In the limit of flux freezing, r„i.o ~~ ► 0, we obtain the 
reduced dispersion relation, 

to 2 + 2TrG ( 7„,ofc(l - A^o 2 ) - ^(Ceff.o + Va,o) = (18) 

The gravitationahy unstable mode corresponds to one of 
the roots of oj 2 < and occurs for /i > 1. The growth 
time for this mode can be written as 

T = - no) 

9 27r[Ga n , (l-Mo~ 2 )(A-AMs)] 1 / 2 ' V ' 

for A > Xms, where 
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Figure 1. Growth time of gravitationally unstable mode (t 9 /<o) as a function of wavelength (A/Lo) f° r models (a) T n ifi/to = 0, (b) 
Tnifl/to = 0.001, (c) T n ifi/to = 0.04, and (d) T n ifi/to = 0.2. Each panel shows time scale curves for models with mass-to-flux ratios 
Ho' = 0.5, 1.0, 1.1, 2.0, and 6.0 (labeled). 
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Figure 2. Left: Minimum growth time of gravitationally unstable mode (r 9 . m /to) as a function of mass-to-fiux ratio (fio)- Right: Length 
scale of most unstable mode (A 9 ,m/£<o) as a function of the mass-to-flux ratio. 



With the inclusion of ambipolar diffusion, the gravi- 
tationally unstable mode still corresponds to one of the 
roots of the full dispersion relation; however since it is a 
cubic function, the expression corresponding to the roots 
cannot be written down as simply as Equation (I19I) and 
therefore are evaluated numerically. 

Figure Q] shows the growth time T g /to as a function 
of the wavelength for four cases of the neutral-ion col- 
lision time Tnifi/to. Each panel shows the dependence 
for several labeled values of /io- Note that these figures 
are presented in normalized form. These different cases 
were chosen to represent various areas within a molecu- 
lar cloud: diffuse regions with high ionization fractions 
(■?m,oAo = 0.001), dense core forming regions with low 
ionization fractions (r n jo/^o — 0.2) and an intermedi- 
ate region that is somewhat between the two extremes 
(i~ni,o/to = 0.04). The flux- frozen case serves as a ref- 
erence point. The value of T n i t o/to — 0.001 was cho- 
sen to correspond to a visual extinction A v = 1 in our 
adopted ionization model for the cloud (see §02 Figured]). 
Specifically, this value corresponds to a column density 
a n = 0.004 g cm" 2 . 

For the flux-frozen case, note the dependence on the 
value of fiQ] there are no unstable, gravitationally collaps- 
ing modes for /io < 1- This is consistent with the dis- 
cussion earlier that for a flux-frozen cloud, only initially 
supercritical clouds can collapse. For the other three 
cases however, the addition of ambipolar diffusion allows 
for unstable, gravitationally collapsing modes to exist for 
regions with /io < 1. For all four plots and all cases of 
Ho shown, note that each curve has a distinct minima. 



This minima represents the shortest growth time (fastest 
growthrate) and corresponding preferred length scale for 
instability. We find the location of these minima for a 
series of values of /to for fixed values of r n j oAo- The 
left panel of Figure [5] shows the minimum growth time 
of the gravitationally unstable mode (Tn,mAo) as a func- 
tion of the critical mass-to- flux ratio (/io). In the limit 
of flux-freezing, the curve shows that for the supercrit- 
ical regime (/io > 1), the growth time for instability is 
short, essentially the dynamical time 2t rs Zq/c s . As the 
mass-to-flux ratio approaches the critical value (/io = 1) 
the growth time scale for instability becomes infinitely 
long. Consistent with the discussion above, this implies 
that in the absence of ambipolar diffusion, it would take 
an infinite amount of time for a region with a critical 
or subcritical mass-to-flux ratio to collapse into a core. 
As the mass-to-flux ratio approaches infinity, this implies 
negligible magnetic support. In this regime, the growth 
time Tg^T is dependent on the critical thermal length scale 
(A T = G<r„,oCeH o) as follows, 



A 



T f,,T 



2^[G ( t„, (A-At)] 1 / 2 



(21) 



The minimum growth time for the unstable mode occurs 
at At,to — 2At- With the addition of ambipolar diffusion, 
the curves show that the subcritical regime has a finite 
growth time scale for instabilities. The magnitude of this 
time scale depends on the degree of ambipolar diffusion, 
i.e., the time scale for neutral- ion collisions. For short 
neutral-ion collision times (r n i,o/^o < 0.2), the time scale 
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Figure 3. Left: Minimum growth time (T 3jm /to) as a function of the dimensionless initial neutral-ion collision time (r n i o/to)' Right: 
Most unstable length scale (A g , m /Lo) as a function of the dimensionless initial neutral-ion collision time. 



for collapse of a subcritical region is 10 - 50 times longer 
than that of a supercritical region; frequent collisions 
between the neutrals and ions will slow the redistribution 
of matter across the field lines. This is the origin of the 
often quoted result that the ambipolar diffusion time is 
« 10 times the free-fall time. As the neutral-ion collision 
time increases (i.e., r„i,o/io > 0.2), the time between 
collisions increases to a point that neutrals rarely /never 
collide with an ion, and therefore have no knowledge of 
the magnetic field present in the molecular cloud. In this 
case, the time scale for collapse of the subcritical regions 
approaches that of thermal collapse (t^.t). 

The right panel of Figure [2] shows the wavelength with 
min imum growth t ime, \„. m /Lo, as a fu nction of /io (see 
also IMortod HMD : ICiolek fc Basull2006ft . There are sev- 
eral features of this graph that are worth noting here. In 
the flux-frozen the mass-to-flux ratio approaches 

the critical value (/io = 1) the wavelength with the min- 
imum growth time approaches infinity. For a non-zero 
neutral-ion collision time, the wavelength with the min- 
imum growth time becomes finite for mass-to-flux ratios 
less than or equal to the critical mass-to-flux ratio. This 
results in a peak in the region of /io ~ 1- Note that the 
location of this peak changes, becoming slightly more 
supercritical as the neutral-ion collision time increases. 

An alternate way to look at this data is to plot the 
minimum growth time and length scales as a function of 
the neutral ion-collision time for a constant mass-to-flux 
ratio. This is shown in Figure |3] The panels show the 
minimum time scale for collapse (T S)TO /to 5 left) and cor- 
responding wavelength (A Sjrn /£o) right) as a function of 



the dimensionless neutral-ion collision time (r^n/to) for 
the five cases of the initial mass-to-flux ratio (/io) dis- 
played in Figure [TJ Focusing on the left panel first, for 
/io < 1.0, the growth time decreases linearly until the 
neutral-ion collision time reaches unity, after which they 
asymptote to the thermal collapse time. Conversely, for 
/io = 2.0, there is a plateau for highly ionized regions 
where the collapse time is longer while for low ionization 
fractions the collapse time is again the thermal collapse 
time. For large enough values of /io, the gravitational 
effects outweigh the magnetic effects and the cloud will 
just collapse on the thermal time scale. Moving to the 
right panel, the general trends shown by all curves indi- 
cate that the wavelength with the minimum growth time 
is a maximum for short neutral-ion collision times and 
a minimum for long neutral-ion collision times. For re- 
gions with high ionization fractions (r n i,oAo ~ 0.1), as 
/ito increases, the wavelength increases from the thermal 
wavelength (At) to a maximum at /io = 1.1, and then de- 
creases back to the thermal wavelength as /io rises further 
(see Figure [J). On the other extreme, for low ionization 
fractions, the variation between the different mass-to-flux 
ratios is almost indistinguishable since the neutrals are 
poorly coupled to the ions and the magnetic field. Note 
however, that for intermediate values of neutral-ion col- 
lision times, T n ifi/to £ 1.0, the preferred wavelength for 
/io = 2.0 becomes larger than those for /io = 1.0 and 
/io = 1.1. This is a result of the complex interplay be- 
tween ambipolar diffusion and field-line dragging in this 
hybrid regime of moderate coupling. 
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3. CLUMPS, CORES AND THE EFFECT OF /j 

Observational studies of star forming regions have 
shown that there are several different extinction thresh- 
olds which define the state of a molecular cloud. 
iRuffie et all ()1998[ ) identify a threshold of A v ~ 3 
whereby above this value clumps within the Rosette 
Molecular cloud are believed to have embedded stars. 

Studies of the Ophi uchus and Pers eus clouds by 

IJohnstone et al.l (|2004H and IKirk et al.l (120061 ). respec- 
tively, suggest that there is a core formation threshold of 
A v ~ 5 and a star formation threshold of A v ~ 7—8 mag- 
nitud es (see also lOnishi et al.l H998; Froc brich fc Rowlei 

The ionization state of a molecular cloud depends on 
its density and consequently the v isual extinction. As 
demonstrated by IRuffle et al.l (|1998D , the ionization pro- 
file of a molecular cloud is one in which the outer lay- 
ers are highly ionized due to photoionization by back- 
ground ultraviolet (UV) starlight, while the inner layers 
are shielded from the UV photons and ionized by cosmic 
rays. This picture assumes a quiescent molecular cloud 
which has no evidence of star formation in its past; ra- 
diation from previous generations of stars can greatly 
complicate this picture. 

In addition to the extinction thresholds, the linear 
analysis discussed above showed that there are correla- 
tions between the mass-to-flux ratio of the region and 
the expected time and length scales for collapse. With 
regards to time scales, subcritical regions evolve more 
slowly than supercritical regions. Conversely, with re- 
gards to length scales, highly supercritical or subcritical 
regions can collapse on smaller length scales than tran- 
scritical regions (i.e., those with fi$ ~ 1). The following 
subsections outline our method and results for our two- 
stage fragmentation model. 

3.1. Model and Calculations 

In their investigation, IRuffle et al.l (|1998| ) , modeled the 
chemical evolution of a collapsing clump using a simpli- 
fied model. Here, we have constructed a model ionization 
profile for the clou d similar to those shown in Figure 1 of 
Ruffl e~et al.l ((1998) . Unlike their ionization profiles how- 
ever, ours takes into account both the UV and cosmic 
ray (CR) ionization regions. The UV region is described 
by a step function of the form 



logx* =logXi,o+0.5(logXi, c -logx i>0 ) x 



1 + tanh 



A — A 



A 



V . d 



,(22) 



where Xi is the ionization fraction, and Xi,0i Xi,c, A VtCr i t 
and A v ^d are the step function parameters, namely the 
maximum ionization, minimum ionization, location and 
width of the step respectively. The cosmic ray regime 
is derived from the well known power-law expression for 
the ion density (rij) as a function of neutral density (n n ), 
namely 

n t = Xn]l 2 , (23) 



where X = lO^cm- 3 / 2 (lElmegreenl 119791: 
lUmebavashi fc Nakano1ll980t lTielensll2005[ ). As shown 
by iCaselli et al.l 1 19981 Figure 11), clear correlations 
between ni and n n are difficult to establish observa- 
tionally. Detailed calculations using grain chemistry 



(|Ciolek fc M ouschovias 1994) show that the exponent 
in the above relation is not a fixed value. However, 
a value of 1/2 is a reasonable theoretical average for 
the relatively low density regime considered here. We 
combine Equation (|23p with vertical hydrostatic equilib- 
rium (Equation [T]). A smooth transition between these 
two regimes is chosen to occur at a visual extinction of 
A v ,cr = 6.365. The full profile is then described by 



0.5(log Xi,c - 

Av-A v ,, 



logXi,o)x 



1 + tanh 



logXi = < 



logXi,o 

( 

log[1.148x 1(T 7 (1 + P cx t) 1/2 x 

f_T\ 1 / 2 ( 2.75 mag ^! 
llO KJ \ A v )\ 



A v < A 



v,CR 



A v > A v ^cr 
(24) 



where the values of the step function parameters are 
logx^o = -4, logXi, c = -7.3203, A v<crit = 4.0 and 
A Vi d = 1-05. Four decimal place accuracy is used for 
some of the profile parameters to ensure a smooth tran- 
sition, rather than for empirical reasons. The midpoint 
of the drop in ionization (A VtCr i t ) has been adopted to 
be closer to A v ~ 5 for comparison to observations in 
the Perseus region (see Section 13.21) . All numbers above 
are to be taken as typical values, and our intention in 
this paper is to outline a general scenario rather than 
performing a detailed parameter search. 

The above profile can also be expressed in terms of 
the column density instea d of the visual extin ction. Fol- 
lowing the prescription of iPineda et al.l (|20100 . the con- 
version from visual extinction to column density can be 
achieved by combining the ratio of H 2 column density to 
color excess (Bohli rTet al.l 119781 ) with the total selective 
extinction (jWhittetl 120031 ) to yield a conversion factor of 
N{H 2 ) = 9.35 x 10 20 (A„/1 mag) cm- 2 . Although this 
conversion is specifically for H2, the abundance ratios of 
other molecules present in the cloud are so much smaller 
(i.e., the abundance ratio of CO to H2 is on the order 
of 10~ 4 ) that they do not contribute significantly to the 
overall number density. Therefore we assume that the H2 
number density is representative of all species. Assum- 
ing a mean molecular weight of 2.33 amu, this translates 
into a mass column density conversion of the form 



a n = 3.638 x 10~ 3 (A,/mag) g cm" 



(25) 



Figure [4] shows the resulting profile for Xi &s a function 
of a n . As shown, the adopted model is one which has 
a steep step function as the ionization fraction drops off 
at A v ~ 3. However, instead of leveling off at the bot- 
tom as the step function would normally, it continues to 
decrease as dictated by the cosmic ray relation. From 
this model ionization structure of the cloud, the neutral- 
ion collision time as a function of column density and 
ionization fraction (Equation [8]) can be derived. 

For each value of r n i^, assuming a constant [Iq, Tg <m 
and X g m can be determined in a similar manner as in 
Figure [3J Figure [5] shows the minimum time scale (top) 
and corresponding length scale (bottom) as a function of 
o~ n for assumed mass-to-flux ratios of 0.5 (dash-dot-dot 
line), 1.1 (dotted line) and 6.0 (dash-dotted line). Also 
plotted for reference is the ionization curve from Figure^] 
(dashed line). Note that in the lower panel, the derived 
values for the (1q — 0-5 and /io = 6.0 curves are almost 
equal and lie nearly on top of one another. As depictedin 
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Figure 4. Model ionization curve for a molecular cloud. Ioniza- 
tion fraction (xi) is plotted against neutral column density (cr n ) 
and corresponding visual extinction (A v ) based on Equation {25}. 
The curve is a composite which captures both the photoionized 
region (A v < 6.365 mag; <j n < 0.0232 g cm -2 ) and the cosmic- 
ray-ionized regions (A v > 6.365 mag; <j n > 0.0232 g cm -2 ). 

the right hand panel of Figure [2j this correspondence of 
values is evident in the locations of the curves at these 
values of the mass-to-flux ratio. 



For a diffuse molecular cloud, (i.e., A v < 3), the mass- 
to-flux ratio can be expected to be either sub- or tran- 
scritical as discussed in Section Q] Figure [5] shows that 
while the length scale for fragmentation of a subcritical 
region is small, the time scale is essentially infinite, there- 
fore direct fragmentation of a diffuse subcritical cloud 
into cores is unlikely. Instead, the more likely scenario is 
the fragmentation of a transcritical cloud. Under this as- 
sumption, the length scale for fragmentation is on the or- 
der of several parsecs and the time scale for collapse is no 
longer effectively infinite (see Figurc[5j dotted line) there- 
fore clumps can form. As the newly formed clump starts 
to contract, the density increases, and consequently the 
visual extinction. This contraction will eventually push 
the clump past a critical threshold A v-cr i t where the 
ionization fraction, and the length and time scales for 
collapse decrease dramatically. When this occurs, high 
density regions within the clump may subfragment and 



10 



10 



10' 



10 



10' 



10' 



10 

io' 

10 



10 



0.00 



f. \ 



Av 
5.50 



11.00 



i\ \\ 

! '-. \ 1 
> ••. \ \ 

\ \ 

\ ■-■ 



H =0.5 

IV* U 
H =6.0 



10 



10 



-6 _ 

10 X 



10 



0.00 0.02 0.04 

a n (g/cm 2 ) 



0.00 



A v 
5.50 



11.00 



!"■ \ 
I- \ 



"T" 



■ •. \ 

! ••• \ 

\ "-v. 

\ i •• 

\ \ '•. 

\ \ ■ 

v \ 



H =0.5 

IV* u 

H„=6.0 



0.00 



0.02 
a n (g/cm 2 ) 



0.04 



10 



10 



10 



io" 6 ^" 



10 



10 



Figure 5. Derived values of minimum growth time scale (r fl , m , 
top) and corresponding length scale A 9jm , bottom) as a function 
of neutral column density. Plots show curves for assumed mass-to- 
flux ratios fiQ = 0.5 (dash-dot-dot line), /io = 1-1 (dotted line) and 
/io = 6.0 (dash-dot line). Also shown is the ionization curve from 
Figure [4] (dashed line). 

start collapsing independently from the clump. The time 
scale for the collapse of these subfragments is dictated 
by the local mass-to-flux ratio. Trans- and supercritical 
subfragments will collapse faster than the parent clump 
while subcritical subfragments will not. 

From these plots, relevant time and length scales cor- 
responding to the collapse of clumps and cores can be 
determined. Table Q] shows these scales for the three 
extinction thresholds and four values of ■ For each ex- 
tinction threshold, the variation of X g . m and r ffjTO across 
the different values of (1q further highlights the trends 
shown in Figure [21 Looking at a specific /io value, the 
length and time scales for collapse become progressively 
smaller for increasing extinction. Of particular interest 
are the comparisons between the length and time scales 
for the formation of clumps versus cores. Assuming a 
clump has fiQ = l.l and A v ~ 3 and a subfragment has an 
A v ~ 5, A c i um p/A corc = 3 22 and T c } um p/T corc — 2.2 14 
for /io > 1.1 in the core. For /io < 1 in the core, 
Tclump/ T core < 1 indicating that a subcritical subfrag- 
ment will not collapse before the parent clump does. 

Note that these time scales represent one e-folding of 
the density. Depending on the initial density, it may take 
several e-folding times for a structure to reach the thresh- 
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Table 1 

Derived time and corresponding length scales for visual extinction thresholds 



Mo 


A v 


= 3 


Ay 


= 5 


Ay 


= 8 


A 9 ,m (PC) 


r g (Myr) 


A 9 ,m (PC) 


r g (Myr) 


A 9 ,m (PC) 


r g (Myr) 


0.5 


0.363 


931 


0.215 


3.07 


0.134 


0.997 


1.1 


4.717 


2.44 


1.520 


1.09 


0.883 


0.528 


2.0 


0.650 


0.44 


0.372 


0.26 


0.224 


0.157 


6.0 


0.368 


0.29 


0.218 


0.17 


0.136 


0.107 



old density that defines it as a core. Taking all of this 
into account, this analysis suggests that a smaller, more 
dense region can form within a larger, more diffuse region 
and can collapse within the time scale of collapse for the 
larger region. The disparity of time and length scales for 
collapse between the clump and core is far greater for a 
model with magnetic fields and ambipolar diffusion than 
it is for a simple hydrodynamic model. 

3.2. Comparison to Observations 

The above analysis reveals several physical conditions 
that are important for the collapse of clumps and sub- 
clumps into cores. First, for regions with transcritical 
mass-to-flux ratios, the preferred length scale decreases 
as the initial column density /visual extinction of the re- 
gion increases. Based on our ionization model (see Fig- 
ure S]), less dense regions (A v < 3 — 5) have larger ioniza- 
tion fractions (smaller neutral-ion collision times) and 
thus yield a larger length scale for collapse while the 
denser regions (A v > 5) have smaller ionization fraction 
(larger neutral-ion collision times) and thus yield much 
smaller collapse length scales. From a time scale point 
of view, it seems favorable for regions within a clump to 
fragment and start to collapse before the overall clump 
structure collapses or fragments. 

Observations of nearby molecular clouds, includ- 
ing (but not limite d to) the Taurus molecular cloud 
(iRebull et all 120101 : iHartmannl J2002t iSchmalzl et al.1 
I2010D Rosette Molecular c loud (iwilliams et all I1995D 
Pipe Nebula (|Frau et al.l l2010tHRoman-Zuniga et all 



20101) and Perseus molecular cloud (jKirk et al.l 120061 : 



Sadavov et al J [20121 ) show clumpy filamentary structure 
both with and without active star formation. In the fol- 
lowing discussion we choose a small sample of clumps 
which show evidence of substructure, and perform in- 
depth case studies comparing our two-stage fragmenta- 
tion model to the observed clump data. The three re- 
gions that we have chosen come from the Perseus Molec- 
ular cloud, the Taurus Molecular cloud, and the Pipe 
Nebula. These case studies are by no means a thorough 
compilation of all possible regions, but rather a represen- 
tation of some interesting cases. This analysis can shed 
some insight into the possible magnetic field strengths 
within these structures. Note that the existence of sub- 
fragments within these chosen areas indicates that the 
clump has evolved to a point in the parameter space of 
Figure [5] that is to the right of Aj.crit, although the ini- 
tial fragmentation event that created the clumps would 
have occurred at lower densities on the left hand side of 

Ay^crit ■ 

3.3. Case Study A: Perseus Bl-E 

The Perseus molecular cloud is located at a distance of 
about 250 pc in the constellation Perseus. It shows clus- 



tered star formation primarily in two regions (NGC 1333 
and IC 348) while elsewhere in the cloud it is rela- 
tively quiescent. A particularly interesting region is 
a ~ 0.1 deg 2 clump ro ughly 0.7 degrees e ast of the 
Bl clump, coined Bl-E (|Sadavov et al.ll2012l ). This re- 
gion is sandwiched between the two star forming regions 
(NGC 1333 to the west, and IC 348 to the east), how- 
ever previous submillimetcr observations of this region 
showed that despite its high extinction (A v > 5), it 
contained neither dense cores nor youn g stellar objects 
(jKirk et all 120061: Uorgensen et al.l 120071 . among others). 
However, recent data from the Herschel Gould Belt sur- 
vey and subsequent observations from the Green Bank 
Telescope (GBT) have r evealed that there i s substruc- 
ture w ithin this regio n dSadavov et al.ll2012f ). In their 
paper, I Sadavov et al.1 (|2012f ) suggest that the delay in 
the formation of cores within this region could be due to 
a strong magnetic field. In this case study, we use our lin- 
ear analysis methods described above to see if theoretical 
models exhibiting the same density as the observations 
result in similar observed length scales. We also com- 
ment on the mass-to-flux ratio that is likely required to 
generate these observed length scales. 

First, the Bl-E clump as a whole was found to have 
a column density N n = (6.3 ± 3.6) x 10 21 cm -2 and a 
radius of 0.46 pc (| Sadavov et al J 120121 ). Taking into ac- 
count the error range on the density, this corresponds 
to a visual extinction range of 2.89 to 10.59 magnitudes, 
with the mid-point value corresponding to a visual ex- 
tinction of ~ 6.75. Table [2] shows the derived length 
scales for collapse to be compared with the diameter of 
the structure, for the low, mid and high values of A v 
for /io = 0.5, 1.1,2.0 and 6.0. From this analysis, we 
see that one scenario stands out as a good fit. The di- 
ameter of the Bl-E clump is ^ 1 pc which matches the 
model with A v = 6.8 and fio — 1.1 (A g m = 1.05 pc). All 
the other models can be ruled out because they result 
in length scales that are much smaller than the diameter 
of the region. The model with A v — 6.8 and /io = 1.1 
corresponds to a growth time of about 0.5 Myr. 

We n ow turn to the substru ctures identified within this 
clump. iSadavov et al.l (|2012l ) list a range of densities for 
each fragment. Table [3] shows the length scales, A ff , m , 
for each fragment. We have calculated A Sj?n for both the 
upper and lower limits of these densities. For clarity, we 
present only the mean value of A ffjTO for each of the four 
mass-to-flux ratios. Looking at the derived values, we 
see that those found for both fio = 0.5 and /io = 6.0 
are consistent with the observed diameter D Q b s for each 
fragment. This suggests that the cores are in either a 
strongly or weakly magnetized region depending on if the 
mass-to-flux ratio is sub- or supercritical respectively. 

The lack of star formation in this region is evidence for 
either a slower fragmentation process or a later fragmen- 
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Table 2 

Derived length scales for Bl-E 











(PC) 




A v (mag) 


N n (10 21 cm" 2 ) 


Mo = 0.5 


Mo = 1.1 


Mo = 2.0 


Mo = 6.0 


2.9 


2.7 


0.373 


4.85 


0.668 


0.378 


6.8 


6.3 


0.160 


1.05 


0.268 


0.162 


10.6 


9.9 


0.102 


0.66 


0.170 


0.103 



Table 3 

Mean derived length scales for fragments within Bl-E 













(PC) 




Object 


N n (10 21 cm" 2 ) 


D obs (pc) 


Mo = 0.5 


MO = 1.1 


MO = 2.0 


Mo = 6.0 


Bl-El 


12.1 - 35.0 


0.066 


0.056 


0.367 


0.093 


0.056 


B1-E2 


13.4 - 30.0 


0.064 


0.054 


0.357 


0.091 


0.055 


B1-E3 


8.9 - 26.8 


0.072 


0.075 


0.493 


0.126 


0.076 


B1-E4 


7.5 - 31.3 


0.056 


0.083 


0.547 


0.139 


0.084 


B1-E5 


6.5 - 24.5 


0.087 


0.098 


0.645 


0.164 


0.100 


B1-E6 


9.8 - 29.1 


0.055 


0.069 


0.453 


0.115 


0.070 


B1-E7 


9.6 - 29.1 


0.055 


0.070 


0.459 


0.117 


0.071 


B1-E8 


6.5 - 25.1 


0.068 


0.098 


0.643 


0.164 


0.099 


B1-E9 


9.1 - 28.2 


0.051 


0.073 


0.481 


0.123 


0.074 



tation event. The former possibility requires the pres- 
ence of a strong magnetic field in the region. Formation 
of subcritical cores within a transcritical clump would re- 
quire fragmentation along the magnetic fi eld lines, which 
is not currently included in our model. ISadavov et al.l 
(2012) note that this region is likely isolated from the 
surrounding regions given it does not fall into the age 
gradient provided by the two closest regions (IC 348 and 
NCG 1333). This lends credence to the possibility that 
this region is younger than the surroundings. However, 
if the clump/core/star formation threshold values hold, 
A v values for the cores of ~ 7-13 mag suggest that they 
should have formed stars. The lack of observed stars 
within this r egion therefore s ugges ts that collapse is be- 
ing slowed. ISadavov et al.l (|2012l ) concluded that this 
delay is due to a strong magnetic field. 

Currently there are no magnetic field measurements 
specifically for the Bl-E region, however measurements 
of the nearby Bl region have shown t he existence of a 
stron g magnetic field (B = 19 — 27/iG) (|Goodman et al.l 
119891: iCrutcher et al.lfl993h . The mass-to- flux ratio of a 
cloud ca n be written in terms of the density and field 
strength ([Chapman et al.ll2011f ) as, 

H = 7.W [l (H 2 )/B t ot, (26) 

where Nn(H2) is the column density in units of 10 21 

cm -2 along a magnetic flux tube and Btot is the to- 
tal magnetic field strength in fiG. Observable quantities 
are the column density of molecular hydrogen, N(H2), 
and the magnetic field component parallel to the line 
of sight (-B||). Projection effects between Nu {H2)/B to t 
and the observed quantity N{H 2 )/B\\ will overestimate 
fi by an average factor of 3 assuming a random orien- 
tation of_tfifi_jna£neH£fiekiwith respect to the line of 
sight (Hcilcs fc Crutcherl 120051 ). By combining this with 
the conversion factor for H 2 number density to visual ex- 
tinction, namely N(H 2 ) — 9.4 x 10 20 A„, the mass-to-flux 
ratio of a region can be calculated in terms of the visual 
extinction and line-of-sight component of the magnetic 



field, 

fi = 2AA v /B ll , (27) 

where A v is in mag and 5m is in /iG. 

Assuming that the magnetic field in Bl-E is similar 
to its neighboring region Bl and applying this equation 
to the values for the cores, we find that the majority 
of the mass-to-flux ratios fall into the transcritical to 
supercritical regimes (0.7-4.42). Some of the lower end 
densities for the cores give significantly subcritical val- 
ues (/x ~ 0.32), however given the range in the observed 
densities for each core, these lower values are likely not 
indicative of the true density. Less extreme mass-to-flux 
ratios (i.e., \i ~ 4 vs [i ~ 6) suggest tha t the c ores may 
be larger than defined in ISadavov et al.l (|2012f ) . Follow- 
ing this scenario, a strong magnetic field combined with 
higher densities results in regions that are transcritical 
which collapse on longer time scales than their super- 
critical counterparts. 

Based on the above analysis, we conclude that this 
region is evolving at a slower pace than the neighbor- 
ing regions and will need to attain greater column den- 
sities than the typical threshold value (A v ~ 8) to 
push the cores into the supercritical regime, due to the 
strong magnetic field p resent. Observa t ions o f the en- 
tire Perseus region by iGoodman et al.1 (|1990l) qualita- 
tively indicate that the magnetic field strengths in NGC 
1333 and IC 348 are weaker than in the Bl and Bl-E re- 
gions. Assuming these observations are fairly complete, 
the presence of a weaker field in the regions which ex- 
hibit star formation and a greater field in those regions 
which do not show signs of star formation are generally 
consistent. 

3.4. Case Study B: Taurus B218 

The Taurus molecular cloud is one of the closest 
star forming region s, at a distance of only 137±10 pc 
(|Torres et al.l 12007). It is known to be host t o at least 
250 young stellar objects (|Rebull et al.ll2010l) that are 
distributed along the filamentary structures within the 
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Table 4 

Derived length scales for B218 









,m 


(PC) 




A v (mag) 


N n (10 21 cm" 2 ) 


Mo = 0.5 


flQ = 1.1 


Mo = 2.0 


Mo = 6.0 


3 


2.80 


0.363 


4.717 


0.650 


0.368 


5 


4.67 


0.215 


1.520 


0.372 


0.218 


10 


9.39 


0.1074 


0.7047 


0.1795 


0.1090 



Table 5 

Derived length scales for fragments within B218 



A g ,m (pc) 

Object N n (10 21 cm" 2 ) D obs (pc) ftp = 0.5 MO = 1.1 fiQ = 2-0 MO = 6.0 



B218-1 


10.4 


0.26 


0.097 


0.633 


0.162 


0.098 


B218-2 


5.9 


0.24 


0.170 


1.116 


0.284 


0.172 


B218-3 


10.9 


0.16 


0.093 


0.609 


0.155 


0.094 


B218-4 


9.4 


0.16 


0.107 


0.705 


0.180 


0.109 


B218-5 


6.5 


0.14 


0.156 


1.020 


0.260 


0.158 


B218-6 


6.5 


0.14 


0.156 


1.020 


0.260 


0.158 



cloud (|Hartmannll2002t ). iSchmalzl etafl (|2010t l present 
a high resolution column density map of the L1495 fil- 
ament. Extinction maps of this region show that the 
lowest visual extinction is on the order of A v = 5 mag 
and numerous dense fragments exist with peak extinc- 
tion values of A v > 15 mag. Specifically, there are 4 
fragments that coincide with fairly large dense core pop- 
ulations. These fragments are four known Barnard Ob- 
jects (B213, B216, B217 and B218) containing 14, 10, 8 
and 6 dense cores respectively. The size of the Barnard 
Objects are ~ 1 pc across while the radii of the cores are 
~ 0.1 pc. For our analysis, we choose to focus on B218. 

B218 is an elongated object with an area of 0.13 
pc 2 . The major axis is 0.82 p c while the minor axis 
is 0.35 pc. Schmal zl et al.l (|2010fl quote a surface density 
N n — 9.39 x 10 21 cm~ 2 for this region, corresponding to 
a visual extinction of about 10 mag. However, this only 
includes regions that have A v > 5 mag. Other observa- 
tions of the region suggest an overall visual extinction o f 
the B218 object to be about 3 mag (G aida et al.l 11984 ). 
For our analysis, we will look at three extinction values, 
namely A v — 3.0, 5.0, and 10. Table |4] shows the de- 
rived length scales A 9iTO for mo = 0.5, 1.1,2.0, and 6.0. 
From this analysis, assuming that the length of the major 
axis corresponds to the collapse length, we see there is 
no exact match between the derived and observed length 
scale for the extinction values used. However, based upon 
the size scales derived for the various mass-to-flux ra- 
tios, we can conclude that the clump is transcritical; the 
length scales for the other mass-to-flux ratios are too 
small under all A v considerations. Specifically, the ob- 
served length scale falls within the parameter space de- 
fined by 5 mag < A v < 10 mag and 0.9 /l«o < 1.1. 
A better match could be found by either increasing the 
visual extinction or decreasing value of the mass-to-flux 
ratio. Given that the mass-to-flux ratio will increase as 
the density increases, we suggest the former possibility; 
the visual extinction of the region is likely between 5 and 
10 mag. Looking back at Table [TJ we see that a region 
with an A v — 8 and /io — 1.1 gives a length scale more 
on par with the major axis of B218. 

From this analysis we can see that by identifying the 



major axis as the collapse length, the mass-to-flux ratio 
of the clump is estimated to be 0.9 < /io < 1.1. Re- 
cent polarization observations of the Taurus region that 
estimate the magnetic field stren gth based on the dis- 
persion of polarization directions ([Chapman et al.ll2011l ) 
show that the entire cloud complex which contains B213, 
B216, B217 and B218 (defined collectively in their paper 
as B213) is significantly subcritical. These two results 
are in direct conflict with each other. The main reason 
for this discrepancy is the lumping together of the four 
Barnard objects into one larger region. Specifically there 
is a discrepancy in th e visual extinction of the regions. 
IChapman etafl poll cite a visual extinction for their 
B213 region of ~ 2 mag while all four of the Barnard ob- 
jects that make up this region are quote d to have visual 
extin ction values on the order of 10 mag (Schmalzl et al. 
120101 ). However, based on the analysis above, we see that 
the actual visual extinction of the B218 region is implied 
to be about 8 mag. Using this value of A v and the mag- 
netic fi eld strengths of B213 measured by Chapma n et al.l 
(|2011f ) in Equation ([27|) results in only somewhat sub- 
critical values (/i ~ 0.7) as compared to those computed 
us ing A; = 2. Also , given that the region defined as B213 
in IChapman et all (j 2 1 If ) covers a much larger area than 
the B218 region, it is plausible that the local magnetic 
field strength in B218 could be greater than the global 
field of their B213 region. Taking errors into account, it 
is possible for these regions to have mass-to-flux ratios 
that are closer to the critical value, as suggested by our 
analysis. 

Turning to the substructures, ISchmalzl et al.l (|2010l ) 
list six high density regions within B218. Table [5] shows 
the derived length scales A ffjJTl for each fragment. In this 
case, there is no clear cut mass-to-flux ratio regime that 
yield a match between X g ^ m and -Dobs for all of the cores. 
Fragments B218 2-4 have diameters that agree well with 
Ag jlrl in the /io — 2.0 regime while fragments 5 and 6 agree 
well with Xg jm in the /io = 0.5 and fj,o — 6.0 regimes. Fi- 
nally, fragment 1 is consistent with a mildly supercritical 
mass-to-flux ratio (1.1 < Mo < 2.0). As with the Bl-E 
region, the lack of observed star formation within B218 
is likely due to the strong magnetic field. Overall, the 
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Table 6 

Derived length scales for B59-09 













(PC) 




A v (mag) 


N n (10 21 cm" 2 ) 


D obs (pc) 


Mo = 0.5 


MO = H 


Mo = 2.0 


Mo = 6.0 


10 


9.35 


0.3 


0.1078 


0.7071 


0.1802 


0.1093 


31.6 


29.6 


0.11 


0.0341 


0.2233 


0.0569 


0.0345 



Table 7 

Derived length scales for select cores within the Pipe Nebula 













(PC) 




Object 


N n (10 21 cm" 2 ) 


Dobs (pc) 


Mo = 0.5 


MO = 1.1 


MO = 2.0 


Mo = 6.0 


14 


13.3 


0.071 


0.076 


0.50 


0.127 


0.077 


40 


11.1 


0.104 


0.091 


0.60 


0.152 


0.092 


48 


6.1 


0.127 


0.164 


1.08 


0.274 


0.166 


109 


47.6 


0.063 


0.021 


0.14 


0.035 


0.021 



analysis of this region yields results that are consistent 
with the results of Section 13.11 

3.5. Case Study C: Pipe Nebula 

Finally, the Pipe Nebula is slightly further away than 
the Taurus molecular cloud, at a distance of 145 pc. Un- 
like the previous two regions, the Pipe Nebula, although 
filamentary, shows little evidence of star for mation; the 
only e vidence being in B59 . Rece nt studies (jFrau et al.1 
[201(1 : 1 oman-Zuniga et all 120101 ) show tha t there ex- 
ists ev olutionary variation across the nebula. iFrau et all 
(2010) propose that the fragmentation in the bowl cor- 
responds to early stages of evolution, the collapsing ma- 
terial in the stem corresponds to an intermediate phase, 
and the star formation within B59 corresponds to the 
latest stage of evolution. This in itself makes the Pipe 
Nebula an interesting observational case on its own, as it 
seems to be providing us with several of the snap shots 
needed to decipher cloud to star evolution in one lo- 
cation. Dust emission maps show that there is a dif- 
ference in the amount of structure within each region 
([Roman-Zuniga et al.ll2010|) that follow s the evolutionary 
gradient proposed bv lFrau et al.1 (|201Q[ ) . The maps of the 
shank/stem regions show that they are much more diffuse 
than the B59 region, suggesting that they are indeed less 
evolved. iRoman-Zuniga et al.l ([2010) find 220 extinction 
peaks within the stem, shank, bowl, and smoke regions of 
the nebula. All of the extinction peaks are within regions 
that have visual extinctions g reater than the threshold 
value quoted by iRuffle et al.l ([1998). In addition, they 
all have radii that are on the order of a tenth of a parsec 
or less, which is in line with the sizes found in our anal- 
ysis. Analysis of the location of peaks within the various 
regions of the cloud show that th ere is a high density 
of pe aks within the bowl and B59 ([Roma n-Zunig a et al.1 
[2010) . 

The previous two case studies have looked at the rela- 
tionship between the length scales for clumps and cores 
in regions with little to no star formation. For this re- 
gion, rather than repeating this type of analysis, we use 
our model to look at two different phenomena. First, we 
will apply our model to B59 as a whole. Se c ond, we ex- 
amine the four cores studied bv lFrau et al.l (|2010). The 
four cores in question are core 14, 40, 48 and 109, which 
come from three distinct regions of the Pipe Nebula; 14 



resides in B59, 40 and 48 reside in the stem and core 109 
is in the bowl. 

B59 is the only region within the Pipe Nebula that 
exhibits active star formati on. It is itself split up into 
several clumps as shown by Roman-Zun iga et al.l (2009, 
Figure 6). Although all clumps show evidence of cores 
within them, only B59-09 harbours visible YSOs. As 
such, B59-09 provides the unique opportunity to look 
at all three levels of evolution w ithin star formation: 
clump , cores and YSOs. Figure 2 of lRoman -Zuhig a et al.l 
(2012) shows the centrally condensed A v profile of 
the B59-09ab cores, centered on the star forming re- 
gion. Assuming a minimum visual extinction thresh- 
old for the region of 10 mag, we find that the B59- 
09 clump has a radius of ~ 30000 AU, which corre- 
sponds to a diameter D b s = - 3 pc. Recent studies 
of B59 by [Roman-Zu niga et al.l ([20121 ) show that the 
star forming core(s) within B59, specifically B59-09a 
and 09b, have a molecular hydrogen surface density of 
Nh 2 = 2.96xl0 22 cm -2 and a diameter D obs = 0.11 pc. 
This surface density corresponds to a visual extinction of 
31.6. 

Table [5] shows the results of our analysis for both 
the extended and local clump regions of B59-09, with 
A v — 10 and 31.6 respectively. Focussing on the ex- 
tended region first and comparing the derived length 
scales to the observed, we conclude that the mass-to- 
flux ratio of the region is somewhere between 1.1 and 
2.0. This means that the star forming clump B59-09 as 
a whole is moving out of the transcritical regime and 
is becoming supercritical. Diving a little deeper into 
this region, the second entry in Table [6] shows that the 
high density region which contains the YSOs, if consid- 
ered separately, is also globally between /io = 1.1 and 
Mo = 2.0. 

Going one step deeper, if the spacing between the 
YSOs (~ 0.06 pc) is assumed to have a near one-to-one 
relationship with the scale of the earlier fragmentation, 
our analysis can shed some light on the probable condi- 
tions in the region at that time. Table [6] shows that for 
both A v values, in order to match the YSO spacing, the 
region is either significantly subcritical or supercritical. 
Evidence of YSOs within the region suggests the latter 
as the former have longer collapse times. 

With that in mind, our analysis implies that in or- 
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der for the YSOs to form with such spacing, the lo- 
cal mass-to-flux ratio in the locations of current YSOs 
must have been larger than the surrounding region to 
promote collapse. Observatio ns show that the magnetic 
field strength in B59 is 17/zG (|Alves et al.H2008f >. We can 
therefore narrow the visual extinction range in which the 
YSO fragmentation occurred. In order for a subregion to 
collapse, the mass-to-flux ratio must be sufficiently differ- 
ent from the background value. If we assume a minimum 
mass-to-flux ratio fj, = 2.0, for B = 17 /j,G, Equation (|27|) 
yields a minimum visual extinction of 14 mag. Therefore, 
for the extinction range, 14 mag < A v < 31.6 mag, the 
observed magnetic field results in a corresponding mass- 
to-flux ratio range for the cores to be 2.0 < /i < 4.5. 

We now consider the four cores from the various re- 
gions in the Pipe Nebula. Table [7] shows the results of 
this analysis. Comparing the derived values of \ g , m with 
-D bs, cores 14, 40 and 48 seem to agree well with those 
in either the subcritical regime or highly supercritical 
regime (/i = 6.0). Core 109 does not agree well with 
any of the three regimes, however it can be concluded 
that it is not extremely sub- or supercritical but rather 
somewhere between /io =1.1 and 2.0 . With regards t o 
the evolutionary gradient proposed bv lFrau et al.l (2010). 
our length scales do show a trend to smaller length scales 
for cores 48, 40 and 14. Core 109 does not seem to fit in 
with this trend at all and given its size compared to core 
40 could be mistaken as a star forming core. Comparison 
between the length scales for B59-09 and core 14 is not 
applicable since the latter is not contained in B59-09. 

The evolutionary gradient between B59 and the bowl 
of the Pipe Nebula is a direct consequence of the different 
magnetic field strengths present in the different regions; 
B = 17,30, and 65 /zG in B5 9, the stem, and the bowl 
respectively (jAlves et a,l.ll2008h . The region with the low- 
est magnetic field has formed stars while the other two 
regions have only formed cores. This is consistent with 
the resu lts of our other two case studies (see Sections 1331 
and I3.4[) . whereby the delay in star formation is due to 
a strong magnetic field which forces the region to attain 
a higher column density before it becomes supercritical 
and can collapse. 

The current state of B59 (i.e., shows evidence of sev- 
eral cores/YSOs yet has not collapsed itself) is sup- 
ported by its apparent transcritical nature. The age 
of the stellar cluster within B59-09 indicates that the 
clump has survived for lon ger than 10 free fall times 
(R oman-Zuniga et all 120121 ) . This again suggests that 
the magnetic field within the region is strong enough 
to delay the overall coll apse of the B59-09 clump 
(|Roma n-Zuniga et al. 2012) while the inner regions have 
attained suffic ient conditions to form stars . Based on the 
current data, iRoma n-Zuhig aet al.l (|2012l ) propose that 
B59-09 is currently contracting toward a denser configu- 
ration. Our model of a transcritical clump is consistent 
with this assertion. 



4. DISCUSSION 

The underlying model and subsequent results for the 
two-stage fragmentation scenario, as outlined in Sec- 
tion 13.11 propose various constraints on the conditions 
required for such a model to work in nature. Namely, the 
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Figure 6. Illustrative plots for two stage fragmentation scenario. 
Panels show the same mass-to-flux ratio curves as Figure(5] Shaded 
regions show the parameter space for the formation of clumps and 
cores respectively. 



clump must form in an initially transcritical diffuse envi- 
ronment while core formation requires a denser supercrit- 
ical region within the parent clump. From the three case 
studies presented above, we have found that the observed 
length scales and corresponding mass-to-flux ratios de- 
termined from our analysis tend to support this scenario. 
In general it seems that parent clumps exhibit mass-to- 
flux ratios that are transcritical (0.9 < /io < 1.1). The 
column densities of these regions place them to the right 
hand side of A v , cr i t (see Figure [5]). This is consistent 
with our model since formation of subfragments requires 
that the clump has already crossed over the critical A v 
threshold. 

Our analysis also shows that the fragments within 
these clumps exhibit length scales that are indicative 
of regions that are either very subcritical or very su- 
percritical. In order to truly delineate between these 
two extremes, high resolution magnetic field data is re- 
quired, however without such data, some conjectures 
about the region can still be drawn. As shown by e.g., 
iCiolek & Mouschoviasl (|1994t l. the mass-to-flux ratio in- 
creases as the density of the region increases. The denser 
subfragments within the clump would be expected to 
have a greater mass-to-flux ratio than the global clump, 
allowing it to collapse faster than the clump itself (see 
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Figure [5]) . It is possible for a minority of sub fragments 
to locally have a subcritical mass-to-flux ratio, however 
these subfragments would likely not be able to collapse 
faster than their surroundings due to the long time scale 
required. Although the lack of magnetic field data for 
most of the regions makes definitively choosing one mass- 
to-flux regime over the other impossible, it can be as- 
sumed that regions with active star formation must be 
within locally supercritical subfragments, while regions 
with little to no star-formation are either too young for 
the subfragments to have collapsed or are largely sub- 
critical. 

Given the variability of the conditions within cloud 
regions, Figure [6] shows an enhanced picture of the two- 
stage fragmentation model. The bottom panel shows the 
fragmentation length scale as a function of column den- 
sity. In the diffuse regions of the clouds, as indicated 
earlier, we expect the mass-to-flux ratio be in the tran- 
scritical regime. The parameter space where we expect 
the fragmentation of the gas into parsec size clumps is 
indicated by the hatched region to the left of A ViCrit . 
As shown in the upper panel, this hatched region cor- 
responds to fragmentation time scales on the order of 
2-10 Myr. As the density of the region increases, as de- 
scribed in Section 13.11 the region will cross over to the 
right hand side of A v>CI ft. The hatched region on this 
side of Au iCr it represents the parameter space for sub-frag- 
mentation within the clump. As shown, a wide range of 
mass-to-flux ratios are considered able to fragment based 
on the decrease of the fragmentation length scale. How- 
ever, the equivalent region in the upper panel of Figure |H] 
shows that although regions within this range of mass- 
to-flux ratios can fragment, only the ones with trans- 
or supercritical mass-to-flux ratios (dotted and dash-dot 
lines respectively) are able to collapse at a rate that is 
faster than that of the transcritical clump itself. 

5. SUMMARY 

We have studied the effect of ambipolar diffusion on 
the formation of stellar clusters using the results of linear 
analysis. By combining the linear analysis with realistic 
ionization profiles for a molecular cloud, our analysis has 
yielded a two-stage fragmentation model for clustered 
star formation that includes the formation of clumps and 
their subsequent subfragmentation. We present several 
interesting concepts that are worth noting. 

1. Linear analysis shows that there are varying length 
and time scales for collapse depending on both 
the mass-to-flux rati o and the neu t ral-ion collision 
time of the region (IMortonl fl99lt ICiolek fc Basul 
120061 this work). Transcritical mass-to-flux ratios 
can have significantly larger fragmentation scales 
compared to super- and subcritical regions (Fig- 
ure |2j right). Subcritical cores have significantly 
longer collapse time scales compared to their trans- 
and supercritical counterparts (Figure [21 left). In- 
creased neutral-ion collision times serve to decrease 
the length and time scales for collapse, particularly 
for the transcritical and subcritical regimes. 

2. A combination of the linear analysis and ionization 
profile shows that molecular cloud conditions are 
favorable for a two-stage fragmentation process for 



the formation of stars. This model suggests that 
an initially diffuse, transcritical (i.e., ^jlq is approx- 
imately unity) cloud can undergo an initial frag- 
mentation on parsec scales and then undergo a sec- 
ond fragmentation event. This second event occurs 
once the density of the clump increases past the 
threshold where the ionization fraction drops dra- 
matically and the length scale for fragmentation 
decreases steeply (Figure [5]). 

3. Comparison with several observed core and star 
forming regions show that the clump sizes are con- 
sistent with their mass-to-flux ratio being transcrit- 
ical, while the core length scales imply they could 
be sub- or supercritical. Lack of star formation 
within the observed regions (excluding B59) could 
suggest subcritical cores. However, analysis of the 
measured magnetic field for Perseus and Taurus 
showed that these regions are at least mildly su- 
percritical. 

4. The threshold values for clump/core/star forma- 
tion as outlined in Section [3] are dependent on 
the region in question and should not be used as 
strict values. As shown in our case studies, re- 
gions with higher intrinsic magnetic fields will re- 
quire a greater column density (higher visual ex- 
tinction) to push the region to a transcritical and 
supercritical mass-to-flux ratio. The Pipe Nebula 
is a perfect example of this, whereby cores in all 
three regions (B59, stem and bowl) have densities 
that exceed the star formation threshold A v ~ 8 
in other clouds, but only the region with the low- 
est magnetic field strength (B59) currently shows 
evidence of young stellar objects. 
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